Application of metasurface description for multilayered metamaterials and an 
alternative theory for metamaterial perfect absorber 
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We analyze single and multilayered metamaterials by modeling each layer as a metasurface with 
effective surface electric and magnetic susceptibility derived through a thin film approximation. 
Employing a transfer matrix method, these metasurfaces can be assembled into multilayered meta- 
materials to realize certain functionalities. We demonstrate numerically that this approach provides 
an alternative interpretation of metamaterial-based perfect absorption, showing that the underly- 
ing mechanism is a modified Fabry-Perot resonance. This method provides a general approach 
applicable for decoupled or weakly coupled multilayered metamaterials. 



Electromagnetic (EM) metamaterials are artificial ma- 
terials that can be engineered to exhibit controlled opti- 
cal properties not found in nature over most of the EM 
spectrum^. They usually consist of multiple identical 
layers of periodically arranged artificial structures, and 
are considered bulk homogeneous media with constitu- 
tive parameters obtained by a retrieval based on effec- 
tive medium theory 3 -*^. Recently, heterogeneous meta- 
materials consisting of two or more distinct layers were 
used to realize functionalities such as perfect absorp- 
tion at THz 5 and infrared frequencies^"—, and EM wave 
tunneling 9 . In that work, each layer of the metamaterial 9 
or all the layers as an entirety 5 7 were considered as a 
homogeneous medium. The effective permittivity and 
permeability, were calculated using an established re- 
trieval procedure 3 -^. However, the metamaterials in these 
systems consist of only one functional layer of artificial 
structures (meta-atoms), which is analogous to a sin- 
gle molecular layer in natural materials. It is challeng- 
ing to define bulk effective permittivity and permeabil- 
ity for such single- "meta- atom" -layer systems, since these 
macroscopic material properties typically result from av- 
eraging fields over many molecular layers. In addition the 
thickness of the effective bulk material in these systems 
is not uniquely defined, which therefore likewise renders 
the effective material properties arbitrar y 1Q i n . Further 
complications arise because metamaterials consisting of 
distinct layers, such as perfect absorbers 5 7 , are inhomo- 
geneous in the wave propagation direction, and cannot 
be strictly considered homogeneous bulk media. 

In this paper, we use an effective medium model that 
treats each layer of the metamaterial as a metasurface 
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with unique effective surface electric and magnetic sus- 
ceptibility, Xse and Xsm- Through a thin film approx- 
imation, we obtain the same equations of Xse and Xsm 
as previous metasurface workii, and also reveal the re- 
lations between surface effective susceptibilities and bulk 
effective material parameters. We then use a transfer 
matrix method to analyze the overall EM properties of 
multilayered metamaterials using the effective material 
parameters (surface susceptibilities) of each layer. We 
find that the overall properties of multilayered metama- 
terials can be determined by their individual layer prop- 
erties, in the absence of inter-layer resonance coupling. 
We also find that such individual layer properties are re- 
sponsible for metamaterial perfect absorbers. This con- 
trasts with previous explanations based on bulk effec- 
tive medium theory 5 -"— . To wit, in previous work, the 
entire metamaterial was considered as a homogeneous 
medium with independently engineered effective permit- 
tivity and permeability to reach the condition e e ff = /i e ff, 
both having large imaginary part resulting in the effec- 
tive refractive index, n = n + in and n ^> n . The 
EM wave thus propagates through the first interface (air- 
metamaterial) without reflection and the strength decays 
rapidly to zero inside the metamaterial before reaching 
the second interface (metamaterial-air). However, our 
results show that the interaction (or the assumed mag- 
netic resonance) between the two metallic layers has a 
negligible effect on the absorption. Instead, the func- 
tional mechanism is the Fabry-Perot interference result- 
ing from the multiple reflections in the cavity bounded 
by two metamaterial layers. Finally, we also find that the 
metamaterial EM tunneling 9 - and the metamaterial anti- 
reflection^ can be explained very well by our approach. 
Our approach is generally applicable for decoupled or 
weakly coupled multilayered metamaterials 13 , where the 
coupling due to evanescent modes is inconsequential. 

We begin with Fig. 1(a), where a single-layer metama- 
terial is considered as a homogeneous thin film with thick- 
ness, d, same as the actual thickness of a single-layered 
of metamaterial structure and approaching zero as com- 
pared to the incident wavelength. Transmission and re- 
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FIG. 1: (a) A schematic of a single layer metamaterial consid- 
ered as a homogeneous thin film and the electric and magnetic 
field across it. E 1 , E r and represent transverse electric field 
of the incident, reflected and transmitted EM wave under nor- 
mal incidence, respectively; Et_, E*+, Ht_ and H*+ are the 
total transverse electric and magnetic fields at each boundary 
of the film; E^ v and H"^ are the average transverse electric 
and magnetic fields inside the film; x se and x sm represent 
the effective surface electric and magnetic susceptibilities; e e ff 
and /i e ff represent the effective permittivity and permeability, 
(b) Schematic of a multilayered metamaterial consisting of N 
layers separated by N-l layers of dielectric spacers. 



flection occurs as a plane EM wave propagates normally 
through the thin film, and generally leads to discontinu- 
ities of the transverse electric and magnetic fields, which 
can be described by the following boundary conditions^: 
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where E £ _ = (1 + R)E\ E, + = TE 1 , H £ _ = (1 - R)H\ 
H t+ = TH\ and H™ are defined in the Fig. 

1(a) caption. The average electric and magnetic fields 
inside the thin film, can be approximately defined as 
= (E £ _ + E £+ )/2 = (1 + R + T)E [ /2 and H* v = 
(H £ _ + H t+ )/2 = (1 - R + T)H { /2 for very thin films, 
i.e., d <C A. e e ff and /i e ff are the effective permittiv- 
ity and permeability of the thin film. The right-hand 
side of equations (1) and (2) contains the bulk mag- 
netic and electric current densities, J m = — kj/i e ffH£ v 
and J e = — icjeeffE^. In the limit d <C A, the thin 
film can also be equivalently considered as a single in- 
terface (metasurface) with surface current density, J se = 
j J e dz = J e d and J sm = J J m dz = J m d, resulting 
from the discontinuity of transverse electric and magnetic 
fields across the thin film, respectively. The surface elec- 
tric and magnetic current densities can be characterized 
by effective surface electric and magnetic susceptibilities, 
J se -iueoXaeE™ and J sm = -iuj^Xsm^-T - We can 
obtain Xse = (e e ff - Xsm = (^eff - where the 
constant 1 results from the permittivity or permeability 
of vacuum when replacing a finite thickness slab by a zero 
thickness surface. Using the previous equations for aver- 
age fields, Xse and Xsm can now be expressed as functions 
of the complex transmission and reflection coefficients T 
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where fco is the wavevector in vacuum. Equations (3) and 
(4) are in consistent with previous work 11 , except for a 
sign reversal for Xsm in Ref. 10, which we believe is a 
misprint. 

Since the transmission and reflection coefficients are in- 
dependent from the thickness, d, the surface susceptibili- 
ties, Xse and Xsm, are well-defined parameters describing 
the properties of single-layered metamaterial in isolation. 
Hence they are distinct from the effective parameters of 
bulk metamaterials, e e ff and /i e ff, which are non-unique 
and depend on the effective metamaterial thickness d. 
Despite this, we also find that the effective permittiv- 
ity and permeability of individual metamaterial layers 
calculated using a usual retrieval procedure^ show good 
consistency with e e ff = Xse/d + 1, and /x eff = Xsm/d + 1 
obtained from equations (3) and (4). The main excep- 
tion is the anti-resonance artifacts obtained from the re- 
trieval procedure and due to periodicity effect are ab- 
sent in the effective surface susceptibilities. This means 
common retrieval procedures may be used to obtain the 
surface susceptibilities of single-layer metamaterials with 
some accuracy. 

Employing a transfer matrix method, we can deter- 
mine the overall transmission and reflection of a decou- 
pled multilayered metamaterial from the effective mate- 
rial parameters derived for each layer. In the following, 
we use a metamaterial perfect absorber in Ref£ as an 
example to demonstrate how to apply this approach and 
reveal its underlying mechanism. 

Figure 2(a) shows a perfect absorber metamaterial 6 
consisting of two layers of metallic structures separated 
by a dielectric spacer. The first metallic structure is an 
array of cross- wire resonators and the second is a metallic 
ground plane. Each can be modeled as a metasurface 
respectively. The whole metamaterial is then considered 
as a three-layered system consisting of two metasurfaces 
separated by a dielectric spacer. 

In the perfect absorber metamaterial, the cross-wire 
structure exhibits electric resonance modes with reso- 
nance frequencies determined by its structural param- 
eters. To obtain the effective material parameters of 
metamaterials, we performed numerical simulations with 
CST Microwave Studio (Computer Simulation Technol- 
ogy GmbH, Darmstadt, Germany), which uses a finite- 
difference time-domain method to determine R and T 
of the metallo-dielectric structures. The unit cell used 
in the simulation for the first layer, MMi, is schemat- 
ically shown as the inset in Fig. 2(b). It consists of 
a gold cross- wire with thickness, d\ = 0.1 /im, width 
w = 0.4 /im, length I = 1.7 jim and period a = 2 /im, 
on the spacer layer with thickness, s = 0.185 /im, and 
dielectric constant, e s = 2.28(1 + 0.04i). Gold is modeled 
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as a Drude metal with a plasma frequency, f p = 2181 
THz, and damping frequency, f T = 6.5 THz 15 . Then the 
effective material parameters of the metamaterial, MMi, 
bounded by vacuum, are calculated using equations (3) 
and (4) with slight modification to handle the substrate 
surrounding the metamaterial structure^. Similarly, the 
gold plate with thickness, d2 = 0.2 /im, was modeled as 
the second layer using a unit cell shown as MM2 in the 
inset of Fig. 2(b). The calculated effective surface elec- 
tric susceptibility of the metamaterial is shown in Fig. 
2(b), where the cross- wires exhibit electric resonances at 
wavelengths of 4.85 fim and 1.68 /im, and the gold plate 
exhibits a plasmonic response with large negative per- 
mittivity over the entire wavelength range from 1.5 to 
8 jim. Importantly, the effective magnetic susceptibility 
(not shown here) for MMi and MM2 was calculated to be 
a constant close to zero over the entire wavelength range. 
Figure 2(c) shows two absorption peaks at the wavelength 
of 1.86 and 6.18 /im, obtained by full EM simulations of 
the entire multilayered metamaterial (solid curve). The 
latter peak corresponds to the absorption reported in 
Ref. 6 . To determine the behavior of the whole absorber 
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FIG. 2: (a) The perfect absorber metamaterial is modeled 
as a stack of three layers, the cross-wire metamaterial, a di- 
electric spacer, and the metallic plate, (b) The real (solids 
curves) and imaginary (dashed curves) parts of the effective 
surface electric susceptibility of the metasurface representing 
cross- wires (red) and metallic plate (blue). The inset MMi 
and MM2 shows the unit cells used in numerical simulations 
to obtain the metasurface parameters, (c) The solid curves 
show the reflectance (blue), i?, and absorptance (red), A, ob- 
tained from direct simulations, while the dashed curves show 
the calculated values, R c (blue), A c (red), using a 3- layer 
metamaterial model. 



we first derive the transfer matrix of the individual lay- 
ers. The transfer matrix of a metasurface, bounded by 
vacuum on each side, can be determined from the relation 
between the transfer matrix and S-parameter matrix, 
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where £21, $12 are forward and backward transmission 
coefficients, and Sn, S22 are reflection coefficients at 
front and back sides of the metasurface, respectively. 
Equation (5) also applies to the dielectric spacer. All 
of the individual transfer matrices are now multiplied to 
obtain the total transfer matrix of the whole metamate- 
rial, M tot =MMM 1 MsM M M 2 . This now constitutes a full 
description of the metamaterial perfect absorber based 
on the effective parameters of the individual metasur- 
faces. Using the relation between the transfer matrix 
and S-parameter matrix again, we can obtain the trans- 
mission^and reflection coefficients of the whole metama- 
terial, T = S21 = l/M\f and R = S n = M\f/M\f.JiS 
shown in Fig. 2(c), the calculated reflectance, R c = \R\ 2 , 
and absorptance, A c = 1 — \R\ 2 — |T| 2 , agree very well 
with the corresponding R and A obtained from a direct 
simulation of the whole metamaterial. Since the transfer 
matrix calculations only take account of the transmis- 
sions and reflections between individual layers, and since 
each layer's properties were determined in isolation, this 
shows that any inter-layer resonance coupling occurring 
between the cross- wire and the metallic plate layers is in- 
consequential. Hence the magnetic response reported in 
previous absorber work 5 7 , is unlikely to have any signif- 
icant functional effect, since it relies on strong coupling 
between two metallic layers in the form of anti-parallel 
resonance currents. To further understand the magnetic 
response, we examined the double-fishnet structure 17 , 
where the magnetic resonance mode exists due to the 
strong coupling between two metallic layers. As we ex- 
pected, the transfer matrix calculation failed to repro- 
duce the magnetic resonance in the direct simulations of 
whole double-fishnet structure since it violates the de- 
coupling or weakly coupling assumption. 

To better understand the mechanism of the perfect ab- 
sorber, we derived the overall reflection coefficient, i?, us- 
ing the mathematical form of a slightly modified Fabry- 
Perot cavity, in terms of the transmission and reflection 
coefficients of metasurfaces: 

(6) 



R = 



R 



1 - R2iR2ze 2{ P 



where, T21, X12, R12 and i?2i, are transmission and re- 
flection coefficients of the metasurfaces, MMi, regarded 
as an interface bounded by semi-infinite media. They 
are generally functions of effective surface electric and 
magnetic susceptibilities of MMi. They can also be ob- 
tained by numerical simulation of MMi using the struc- 
ture shown in Fig. 2(b). R23 = — 1 is the reflection co- 
efficients from the gold ground plane, MM2; /3 = n s kd s 
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is the propagating phase term in the spacer; and a = 
^21^12 — R12R2I' At the perfect absorbing wavelength, 
the reflection coefficient, R = 0, which requires the fol- 
lowing conditions: 

\Ri2\ = \a\ (7) 
<j>{R 12 ) - <!>{ol) - 2/3 = 2ra7r, \m\ = 0, 1, 2, ... (8) 

To understand these conditions, similar to recent work 





Wavelength (mm) 

FIG. 3: (a) The amplitude of R12 (blue) and a (red), (b) The 
phase terms in equation ©, = <j>{R\2) — — 2/3 (red), 
the phase of R12 (blue) and a (green), and the propagation 
phase 2/3 (black) are shown respectively. 

on anti-reflection metamaterials 12 , we calculated the am- 
plitude and phase terms shown in equations ([7j) and (JSj) 
using the effective material parameters of the metasur- 
faces. As shown in Fig. 3(a), in the strongly absorb- 
ing regions (shaded) centered at the wavelengths of 6.18 
and 1.86 /xra, the amplitudes of R12 and a are almost 
equal, roughly satisfying the amplitude condition. In 
Fig. 3(b), the phase term, 6 crosses zero and 2tt at 
wavelengths of 6.18 and 1.86 /im, respectively, indicat- 
ing the phase condition in equation (jSJ is perfectly ful- 
filled at the absorption peaks. Several other absorption 
peaks (not shown here) can be observed at shorter wave- 
lengths when reaches 47r, 6tt etc. Figure 3 also shows 



that equations ([7j) and ([8]) are not simultaneously sat- 
isfied at the same wavelengths, which explains why the 
reflectance \R\ 2 (blue dashed curve in Fig. 2(b)) does 
not reach zero. Equations (|6]l8|) and Fig. 3 indicate that 
the nature of the absorber is Fabry- Perot- like resonance 
modes resulting from multiple wave reflections between 
metasurfaces MMi and MM2. The metasurface MMi 
and ground plane MM2 form a Fabry-Perot cavity filled 
with a lossy dielectric spacer. Strong absorption occurs 
as the EM wave propagates through the lossy dielectric 
spacer multiple times. The metasurface MMi provides 
the proper surface susceptibility to fulfill the conditions 
presented in equations ([7j) and J5J). Practically, the un- 
derstanding of this mechanism helps us to improve the 
metamaterial designs. For instance, to achieve a perfect 
absorption, we need to optimize the design of the meta- 
material, MMi, to reach the conditions in equation 
and (|8]) simultaneously. Adjusting the slope of \R12\j \ a \ 
and curves shown in Fig. 3 can also optimize the ab- 
sorbing bandwidth of the whole metamaterial. 

In conclusion, we have proposed an effective medium 
model for individual metasurfaces. In our model, each 
layer of the multilayered metamaterial is considered as 
a metasurface. These metasurfaces may be combined 
to determine the properties of multilayered metamate- 
rials using the transfer matrix method. This alternate 
interpretation resolves the problems of defining a mul- 
tilayered metamaterial as a single- "meta- atom" -layered 
bulk medium. This method provides a general approach 
applicable for any decoupled or weakly coupled multilay- 
ered metamaterials. We applied this method to the re- 
cently demonstrated perfect absorber metamaterials and 
identified the underlying mechanism as Fabry-Perot type 
resonance modes in contrast to the previously reported 
mechanism of independent engineering of the bulk effec- 
tive permittivity and permeability. We have also found 
that this model accurately reproduces the previously re- 
ported EM wave tunneling effects 9 . 
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